Quantum phases of atomic boson-fermion mixtures in optical lattices 
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The zero-temperature phase diagram of a binary mixture of bosonic and fermionic atoms in an 
one-dimensional optical lattice is studied in the framework of the Bose-Fermi-Hubbard model. By 
exact numerical solution of the associated eigenvalue problems, ground state observables and the 
response to an external phase twist are evaluated. The stiffnesses under phase variations provide 
measures for the boson superfluid fraction and the fermionic Drude weight. Several distinct quantum 
phases are identified as function of the strength of the repulsive boson-boson and the boson-fermion 
interaction. Besides the bosonic Mott-insulator phase, two other insulating phases are found, where 
both the bosonic superfluid fraction and the fermionic Drude weight vanish simultaneously. One 
of these double-insulator phases exhibits a crystalline diagonal long-range order, while the other is 
characterized by spatial separation of the two species. 



Following the seminal experiments on the superfluid to 
Mott-insulator transition in ultracold atomic Bose gases 
in optical lattice potentials 0, 0] , a new frontier in the 
field of degenerate, strongly correlated quantum gases has 
emerged. It has recently become possible to prepare de- 
generate mixtures of a bosonic and a fermionic atomic 
species in an optical lattice 0. These degenerate gases 
will facilitate the experimental study of quantum phase 
transitions in systems of mixed quantum statistics, which 
are extremely hard to access in the solid-state context. 
The superior experimental control available in ultracold 
atomic gases experiments promises detailed insights into 
the rich physics of quantum phase transitions. In this 
paper we explore the phase diagram of a binary boson- 
fermion mixture in an optical lattice and identify exper- 
imental signatures of the different phases. 

A theoretical framework for degenerate boson-fermion 
mixtures in optical lattices can be constructed by a 
straight-forward generalization of the single-band Bose- 
Hubbard model, used successfully to describe the super- 
fluid to Mott-insulator phase transition in purely bosonic 
systems Q. A complete single-particle basis within 
the lowest Bloch band is given by the Wannier func- 
tions for the vibrational ground state of each lattice 
well. We define creation operators &\ (c\) which create 
a boson (fermion) in the localized Wannier state at site 
i = 1, Using this basis we can translate the many- 
body Hamiltonian into a second quantized form. This 
leads to the Bose-Fermi-Hubbard (BFH) Hamiltonian for 
a one-dimensional translationally invariant lattice 0,0,0 



where nf = ajaj and hf = cjfii are the occupation num- 
ber operators for bosons and fermions, respectively. The 
first line describes the tunneling between adjacent lattice 
sites. The tunneling amplitudes J B and Jp are connected 
to off-diagonal matrix elements of the kinetic energy and 
the lattice potential p| . The second line of Eq. JQ) con- 
tains the on-site boson-boson and boson-fermion inter- 
actions with interaction strengths Vbb and Vbf, respec- 
tively. Within the single-band approximation a fermion- 
fermion on-site interaction does not arise as the Pauli 
principle prevents two identical fermions from occupying 
the same site. 



We exactly solve the eigenvalue problem of the BFH 
Hamiltonian within a complete basis of Fock states of the 
form \{nf,...,rf} a )® |{nf , nj}p). By {nf,... ,nf} a 
we denote a set of boson occupation numbers for the in- 
dividual sites. The index a = 1,...,Db labels all pos- 
sible compositions of occupation numbers which yield 
J2i=i n f — Ab. For the set of fermion occupation num- 
bers {nf , ...,nf}p only those compositions with or 1 
fermions at a given site are allowed by the Pauli prin- 
ciple. The dimensions of the boson and fermion space 
are given by D B = (N B + I - 1)\/(N B \(I - 1)!) and 
D F = I\/(N F !(I - A F )!), respectively. For a lattice with 
1 = 8 sites, Ab = 8 bosons, and Ab = 4 fermions the ba- 
sis dimensions are D B = 6435 and Dp = 70. Clearly, the 
dimension of the fermion space is dramatically reduced 
due to the Pauli principle. 
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Within the basis of Fock-states we construct the ma- 
trix representation of the BFH Hamiltonian assuming pe- 
riodic boundary conditions. The lowest eigenvalues and 
eigenvectors of the Hamilton matrix are computed using 
a refined Lanczos algorithm. The eigenvectors provide 
the coefficients for the representation of the eigenstates 



2 



\^ v ) in the Fock-state basis 
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|{n?,...,n?} a )® |{n?,...,nj}/J> • 

From this representation of the ground state we can ex- 
tract various observables such as mean occupation num- 
bers and number fluctuations for both species. It also 
provides access to more elaborate quantities such as the 
one- or two-body density matrix and the static structure 
factor, which we will discuss below. 

The dynamical response of the system to external 
perturbations is a crucial quantity in distinguishing be- 
tween different quantum phases, for example, between 
the bosonic Mott-insulator and the superfluid phase. By 
imposing a linear variation of the phase onto the many- 
body wavefunction one can probe the mobility of the 
atoms in the lattice. On the single particle level, a spa- 
tial variation of the phase of the wavefunction is associ- 
ated with a velocity field vix) — ~V0(x), and a linear 
phase variation thus corresponds to a constant velocity 
across the lattice. Roughly speaking, those particles free 
to move will respond to the phase twist and exhibit a 
homogeneous flow. This flow is then associated with an 
additional kinetic energy and hence an increase of the 
total energy. This provides a measure for the density of 
mobile particles in the lattice 0, |t| . 

This phase variation can be realized by imposing 
twisted boundary conditions for the many-body wave- 
function. In the case of the BFH model it is, however, 
more convenient to map the phase twist by means of a 
unitary transformation onto the Hamiltonian This 
leads to a "twisted" Hamiltonian which contains Peierls 
phase factors in the hopping terms: 
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In order to probe the response of the bosonic and the 
fcrmionic species separately, we have introduced different 
total twist angles Ob and Of, respectively. 

For an isolated boson twist (Ob > 0, Of = 0) the dif- 
ference in the ground state energies of the boson-twisted 
Hamiltonian He B .o and the initial Hamiltonian Ho,o pro- 
vides direct information on the superfluid density of the 
bosonic component. This is because the superfluid com- 
ponent responds to the imposed phase gradient produc- 
ing flow and the associated energy change. The bosonic 
phase stiffness 
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is, therefore, a measure for the boson superfluid frac- 
tion jl, 0, E3> El • We note that neither the stability of 
the superflow nor the reduction of the flow by the lat- 
tice itself are taken into account by this quantity. The 
superfluid weight / S B vanishes for insulating phases such 
as the Mott-insulator phase and goes to 1 for a perfect 
superfluid. 

The energy change resulting from an isolated fermion 
twist (Ob = 0, Of > 0) provides a measure for the mobil- 
ity of the fermionic atoms. The fcrmionic phase stiffness 
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is equivalent to the well-known Drude weight for charged 
fermions in a lattice 0, El> The vanishing of 

the Drude weight is an indicator of an insulating state, 
whereas a finite value indicates non- vanishing conductiv- 
ity. 

These two phase stiffnesses are an important tool to 
characterize the different quantum phases of the boson- 
fcrmion mixture in a lattice. In the following we map out 
the phase diagrams for different boson and fermion num- 
bers, 2Vb and N F , as function of the strength of the repul- 
sive boson-boson and the boson-fermion interaction. For 
simplicity we examine the case of equal tunneling rates 
for the two species, i.e. J = Jb = Jf, and use J as energy 
unit. The only remaining parameters of BFH Hamil- 
tonian are then the ratios between the two interaction 
strengths and the tunneling rate, Vbb/ J an d Vbf / J- For 
each point in the VeB-VeF-plane we solve the eigenvalue 
problem of the BFH Hamiltonian JIJ and the two twisted 

Hamiltonians He B ,o and Ho^Of (with Ob, Of = 0.1) in 
order to compute the phase stiffnesses and several simple 
ground state observables. 

First we consider a lattice with 1 = 8 sites and equal 
filling factors N B /I =1/2 and N F /I =1/2 for the 
bosonic and the fermionic species, respectively. Although 
the system is rather small we could confirm that finite 
size effects are insignificant for the observables discussed 
in the following Figure n shows contour plots of the 
boson stiffness / S B , i.e., the superfluid fraction within the 
bosonic species, and the fermion stiffness fj, i.e., the 
Drude weight for the fermionic component. In addition 
the largest coefficient C 2 lax of the expansion of the 
ground state is given. Large values of C^ ax indicate that 
the ground state is dominated by one or a few particular 
Fock states. 

Figures da) and (b) reveal that in the large param- 
eter region for Vbb/ J > 5 and Vbf/ J > 15 the boson 
and the fermion phase stiffnesses vanish. The bosonic 
and the fermionic components then behave like insula- 
tors, i.e., the superfluid fraction and the Drude weight 
are both zero. This double-insulator phase is gener- 
ated by the combined action of the boson-boson and 
the boson-fermion repulsion. Any configuration which 
has two bosons or a boson and a fermion at the same 
site is energetically unfavorable, and the ground state is 
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therefore dominated by Fock states with either one bo- 
son or one fcrmion per site. Only if one of the interaction 
strengths is small can a finite stiffnesses for bosons and 
fermions emerge. 

Closer inspection of the largest coefficient of the ex- 
pansion J3J, displayed in Fig. dc), reveals that within 
the double-insulating region there are two distinct areas 
where the largest coefficient reaches C^ lax rj O.f . Hence 
there is a class of Fock states which provide the domi- 
nant contribution to the ground state. In region (A) for 
values 2Vbb ^ Vbf the dominant Fock states exhibits an 
alternating occupation of the sites with one boson or one 
fermion. In region (B) the dominant Fock states con- 
tain a continuous block of bosons followed by a block 
of fermions. Note that although these states have the 
largest weight, other Fock states do also contribute to 
the ground state. Most important are those which can 
be generated from the dominant Fock state by exchange 
of two lattice sites or particle-hole excitations. 

The reason for the existence of and the cross-over be- 
tween these two different intrinsic structures is quite sub- 
tle. All the interaction matrix elements for the dominant 
Fock states are irrelevant, because all these states have 
exactly one particle (boson or fermion) at each site. The 
existence of the two different structures is triggered by 
states with double occupancies that enter into the ground 
state with much smaller weight but nevertheless couple 
to the dominant states via the tunneling term. This com- 
plex correlation dictates the cross-over between alternat- 
ing and separated species. 

A good observable to distinguish the different intrinsic 
configurations is the static structure factor 0. For the 
bosonic component it is given by 

S B (ka) = e^-^olnfhf |* ) , (6) 

B i,j=l 

where a is the lattice spacing. It measures the density- 
density correlations in the system and thus provides in- 
formation on the presence or otherwise of diagonal long- 
range order. 
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FIG. 2: Static structure factor S B {ka) for I = 8 and 7V B = 
Nf — 4 as function of ka for the ground state with alternating 
occupation at Vbb/J = 40, Vbf/ J = 15 (solid line) and with 
separated species for Vbb/J = 10, Vbf/ J = 40 (dashed line). 

Figure [21 shows an example for the structure factor 
within the region of alternating occupation (full line) and 
the region of phase separation (dashed line). Apart from 
the trivial peaks at integer multiples of 2n the structure 
factor differs substantially within the two phases. For 
the alternating occupation phase, S B (ka) exhibits an en- 
hancement around ka = it which signals an increased 
probability for a boson to be found at every other site. 
This is the signature of diagonal long-range order as it is 
found in a crystal. In a two-dimensional lattice the cor- 
responding checkerboard configuration, which exhibits 
crystalline diagonal long-range order, was identified re- 
cently 0- 

Within the separated phase the structure factor shows 
a rather different behavior (dashed line in Fig. it 
is suppressed around ka = tt and enhanced at ka w 
2tt/I. Due to this difference in the structure factors it 
might be possible to distinguish the two types of double- 
insulator phases experimentally, e.g., through selective 
Bragg diffraction of light off the trapped atoms. We note 
that the interference pattern produced by the atoms af- 
ter release from the lattice, along with ballistic expan- 
sion, will not distinguish the phases as their coherence 
properties, i.e., the off-diagonal long-range order, are too 
similar. 
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As a second example we consider a lattice with 1 = 8 
sites, a boson filling factor Nb/I = 1, and a fermion fill- 
ing factor Np/I = 1/2. Contour plots of the bosonic and 
fermionic stiffness and the largest coefficient as function 
of Vbb/J and Vbf/J are shown in Fig. |2| 

The qualitative difference compared to the half-filling 
case is the appearance of a bosonic Mott-insulator phase 
(M) for Vbb > Vbf- In this region the boson stiffness, i.e. 
the superfluid fraction, vanishes whereas the fermionic 
stiffness remains at its noninteracting value as depicted 
in Figs. Ota) and (b). Under the influence of the strong 
boson-boson repulsion the bosons form a Mott insulating 
layer. The fermions are not affected, they experience a 
homogeneous background of bosons on which they can 
move freely. 

The situation changes if the boson-fermion interac- 
tion becomes comparable in strength to the boson-boson 
interaction. Through the boson-fermion repulsion the 
fermions are able to break up the bosonic Mott-insulator 
and partially restore bosonic superfluidity in the region 
Vbf ~ Vbb • Further increase of the boson-fermion repul- 
sion leads to the simultaneous reduction of both phase 
stiffnesses and to the formation of a double-insulator 
regime as observed in the case of half-filling. The con- 
tour plot for the largest coefficient C^ ax depicted in Fig. 
|^c) reveals that both of the previously observed intrin- 
sic configurations appear again: (A) an alternating oc- 
cupation of the sites associated with crystalline diago- 



nal long-range order and (B) continuous blocks of bosons 
and fermions which resemble a spatial separation of the 
two components. Due to the appearance of the bosonic 
Mott-insulator both phases have been shifted to the re- 
gion Vbf > Vbb- 

In summary, we have mapped out the phase diagrams 
of a binary mixture of bosonic and fermionic atoms in 
a one-dimensional lattice. By exact numerical solution 
of the eigenvalue problem for the Bose-Fermi-Hubbard 
Hamiltonian we are able to compute ground states as 
well as the dynamical response of the system. The stiff- 
nesses of the bosonic and fermionic component under 
phase twists provides an important measure for the su- 
perfluid properties of the bosonic component (superfluid 
fraction) and the conductivity of the fermionic species 
(Drudc weight). Using these quantities we have identified 
several distinct quantum phases, e.g., a bosonic Mott- 
insulator phase and two double-insulator phases which 
exhibit different degrees of diagonal long-range order — 
alternating occupation and component separation. This 
shows that atomic boson-fermion mixtures in optical lat- 
tices constitute a versatile tool for the study of funda- 
mental quantum phase transitions. 
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